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Abstract 


We present systematic first-principles studies of occupancy, diffusion, and 
segregation of He atoms, as well as underlying evolution mechanism of He 
clusters in Ta, Re and Os precipitates to help to understand the effect of al- 
loying solutes on He retention and bubble formation. It is found that single 
He atom prefers tetrahedral (T) site in Ta and basal octahedral (BO) site 
in Re and Os. Single He atom prefers the diffusion path from one T site to 
the closest T site in Ta, and the diffusion path from one BO site to the next 
BO site passing an octahedral site along (0001) direction in Re and Os. The 
diffusion of single He atom in Ta, Re and Os is relatively difficult than that 
in W. He atoms are energetically favorable to segregate at interstitial sites 
forming He, clusters in Ta, Re and Os, and form close-packed He monolayer 
structure between (110) planes in Ta. The binding strength of He,, cluster 
in Ta, Re and Os is relatively small than that in W (W>Os>Re>Ta). Com- 
pared with interstitial sites, He atoms are more energetically favorable to 
aggregate in vacancy cluster forming Vac,,He,. The binding strength order 
of Vac,,He,, cluster in the systems obeys W~Os>Re>Ta. We thus expect 
that the presence of Ta and Re may contribute to reduce He retention and 
bubble formation in W. 
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1. Introduction 


Confinement thermonuclear fusion energy, a kind of clean and sustainable 
energy source, is expected to provide sufficient energy to satisfy the mounting 
demand of mankind. The International Thermonuclear Experimental Reac- 
tor (ITER) is thereby constructed to realize self-sustaining combustion and 
demonstrate the scientific and technological feasibility of fusion energy [1]. 
Tungsten (W) is chosen as diverter plasma facing materials (PF Ms) in ITER 
and considered a promising candidate for an application in future nuclear fu- 
sion reactors because of the good properties of high resistance to erosion and 
sputtering, and stability against neutron damage. However, the irradiation 
of D, T, and helium (He) ions, and high-energy neutrons (up to 14 MeV) will 
result in the degradation of the mechanical properties and structural strength 
of W by producing H/He bubbles and transmutation element precipitates in 
the material. 

Under neutron irradiation in a fusion reactor, W undergoes transmutation 
to its near-neighbors in the periodic table. The burn-up calculations showed 
that pure W is primarily transmuted to tantalum (Ta), rhenium (Re) and 
osmium (Os) isotopes under high-energy neutron irradiation. By calcula- 
tions, Gilbert and Sublet predicted that the concentrations of transmutation 
elements Ta, Re and Os reach up to 0.81 at.%, 3.8 at.% and 1.4 at.% after 
five years of fusion operation [2], respectively. The concentrations of trans- 
mutation products are unquestionable at significant levels. In contrast, the 
level of transmutation predicted for ITER conditions is lower by orders of 
magnitude [2]. On the other hand, Re is often artificially introduced as alloy- 
ing elements for improving the mechanical properties of W for fusion reactor 
applications, and Ta and Re elements are added in W for reducing the D re- 
tention and suppressing bubble formation. Typically, between 5 and 26 at.% 
Re are added to increase the low temperature ductility and to improve the 
high temperature strength and plasticity, the creep strength and the recrys- 
tallization temperature of W [3-7]. Many experimental works are performed 
to access the influence of Ta and Re on H retention and bubble formation in 
tungsten [8-11]. Zayachuk et al. found that the amount and size of bubbles 
in W-1 at.% Ta and W-5 at.% Ta exposed to D are significantly smaller than 


that in pure tungsten [8-10]. Meanwhile, Alimov et al. found that the size 
and density of bubbles in W-5 at.% Re is much smaller than that in pure 
tungsten exposed to D ions at 523 K [11]. 

Neutron irradiation also induces structure modification of W alloys (W- 
Ta, W-Re, and W-Os) and solute-rich precipitates formation although Ta, Re 
and Os are far below the solubility limit. It was reported that the solute-rich 
precipitates with y phase structure (Re;W) appear in neutron irradiated 
W-based binary alloys containing 5, 11, and 25 wt.%Re [12]. It was also 
found that the x and o phase (ReW) precipitates are formed in W-26 at.% 
Re alloys [13-16]. Similarly, many needle-like precipitates and black dots 
appeared in W-3 at.%Os alloy irradiated to 1.54 dpa [17], and plate or needle- 
like precipitates in {110} planes were observed as well in W-5 at.%Re-3 
at.%Os and W-25Re-30Os dilute ternary alloys [18, 19]. Therefore, there may 
be plenty of solute precipitates in W in a fusion environment. Meanwhile, 
He ions are ubiquitous in W exposed to He ions in a fusion environment, 
and consequently interact with solute precipitates, which will also reversely 
influence He bubble formation and He retention in W. Thus, it is crucial to 
study the behaviors of He in solute precipitates. 

Solute precipitates in W have the common character, i.e., the local con- 
centrations of solutes like Ta, Re and Os are very high. According to the 
phase diagram, W-Re alloy is initially body-centered cubic (bcc) structure 
and transformed to bec+o, o, x and ultimately close-packed hexagonal (hcp) 
phase with the increase of solute concentration [20]. And it is found in our 
previous work that Re and Os are energetically favorable to accumulate at 
vacancies, providing local high-concentration of Re and Os in W [21]. Thus, 
the precipitates with locally high-concentration of solutes are assumed to be 
bec Ta, hep Re and hcp Os. The behaviors of He in solute precipitates re- 
garding the influence of Ta, Re and Os on He bubble formation at atomic 
level are extremely difficult to be directly determined by experiments. While 
theoretical simulations are very effective ways to obtain the behaviors of He 
in Ta, Re and Os. First-principles theory is a kind of accurate and reliable 
method which is widely used to study the energetics of point defects and the 
interaction between He and point defects [22, 23]. In the present work, first- 


principles calculations are systematically performed to study the clustering 


of He atoms at interstitial sites and the interaction of He with vacancy, as 
well as the diffusion behaviors of He in bcc Ta and W, as well as hcp Re 
and Os. It is found that He atoms are energetically favorable to segregate at 
both interstitial sites and vacancies in the systems, and He atoms are more 
preferable to nucleate at vacancies. He atoms prefer to form close-packed 
monolayer structure between (110) planes in bec Ta. The binding strength 
of He, clusters in Ta, Re and Os is relatively small than that in W, and the 
strength of Vac,,He, in Os is close to that in W but larger than that in Ta 
and Re. 


2. Computational details 


All the calculations performed in this work are based on density functional 
theory as implemented in the Vienna Ab-initio Simulation Package (VASP) 
(24, 25]. The projector augment approximation wave potentials (PAW) and 
the generalized gradient approximation (GGA) with the exchange-correction 
functional proposed by Perdew and Wang are employed [26]. The wave 
functions are expanded in the plane-wave basis set using an energy cutoff of 
500 eV for Ta, W, Re, and Os systems. The integrations over the Brillouin 
zone are performed considering a mesh of 3x3x3 k for the calculations in 
the 128-atom 4x4~x4 supercells of Ta, W, Re, and Os. The positions of 
the atoms and the volume of the supercell are fully optimized to obtain the 
energy for each system. When the forces on each atom converge to a value 
below 0.01 eV/A, the structural optimization process will be terminated. 
The point-defect formation energy is calculated according to the following 
formula [27]: 

Ey(X) = E(X) — E(noX) — X` nis, (1) 


where E(X) and E(noX) are the total energies of the systems with and 
without point-defect X, respectively. u; is the chemical potential of species 
i, and n; is the number of atoms of species i added to (n; > 0) or removed 
from (n; < 0) the supercell to form the point-defect. In the calculations, 
the energy of a single atom in a 128-atom supercell is taken as the chemical 
potential for Ta, W, Re or Os, and the chemical potential of helium is zero 


since He is a closed-shell atom. The incremental binding energies of He, 


4 


clusters are then given by: 
E,( Hen) = E(Hen-1) — E(Hen) — (E(P) — Ey(He)), (2) 


where E,(He) is the total energy of the system containing a He atom at the 
energy-lowest interstitial site, and Æ(Hen) is the energy of the system with 
an interstitial He, cluster, and E(P) is the energy of perfect system. n is 
the number of He atoms. The incremental binding energy of VacmHen is 
expressed by: 


E(VacmHen) = [E(VacmHen-1) + Ey(He)| — [E(VacmHen) + E(P)], (3) 


where E(V ac, He,,) is the total energies of the systems containing a vacancy- 
helium cluster Vac,,He,. m=0 or 1. Since He is a kind of light element, 
zero-point-energy (ZPE) correction may influence the stability of He in all 
the systems. ZPE of He atoms are calculated by summing up the vibrational 


energies of the normal modes by the expression: 
1 
ZPE => 2 hy;, (4) 


where fA and v; are Plank’s constant and the normal vibration frequencies. 
The direct calculations regarding the vibrational energies of the systems in 
the presence of He-involved defects is computationally expensive. Besides, 
the masses of the considered metal atoms are much larger than that of He, 
which results in that the vibrational energies of metal atoms are much smaller 
than that of He atoms. Therefore Einstein approximation is used to calculate 
ZPE of He atoms in all metals, and the metal atoms are fixed and ZPE is 


calculated only from normal vibrational modes of He atoms. 


3. Results 


3.1. Solution of interstitial He in Ta, W, Re, and Os 
The most-stable states of both W and Ta are bcc structures, and both 


Re and Os are hep structures. In bcc structure, there are two kinds of 
high-symmetry interstitial sites, i.e., tetrahedral (T) site and octahedral (O) 


site. In contrast, there are as many as six possible interstitial sites in hcp 
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T BT BC BO O C 
Figure 1: (Colour online) The configurations of high-symmetry interstitial sites in hcp 
structure: tetrahedral site (T), basal tetrahedral (BT) site, basal plane (BC), basal oc- 
tahedral (BO) site, octahedral (O) site and the midway between two nearest neighbor 


atoms (C) out of the basal plane. The gray and green balls are Re (or Os) and He atoms, 
respectively. 


structure, that are, T site, basal tetrahedral (BT) site, midway between two 
nearest neighbors in the basal plane (BC), basal octahedral (BO) site, O site, 
and midway between two nearest neighbor atoms out of basal plane (C). All 
the notations of interstitial sites are the same as those in Refs [28, 29], and 
these initial interstitial sites in hcp structure are shown in Figure 1. The 
formation energies of He at the T and O in bcc W and Ta, and He at T, O, 
C, BT, BC, and BO sites in hcp Re and Os are determined. The formation 
energies are calculated with and without ZPE corrections, and the results are 
summarized in Table 1. According to our calculations, He is energetically 
more favorable to occupy T site rather than O site by 0.21 eV and 0.31 eV in 
W and Ta, respectively. The site-preference of He in Ta and W agrees well 
with previous results [30]. The formation energies with ZPE correction for He 
at T and O sites also agree well with previous values [31]. In hcp structure, 
it is unstable for He to occupy BC and BT sites in both Re and Os. When 
He atoms are put at BC and BT, they will move to the closest BO and T 
sites, respectively. Except for the calculations of He atoms at the interstitial 
sites, the substitutional sites are also considered in all the four metals. In 
calculating the formation energy of He at substitutional site, a vacancy is 
assumed to presented in the systems. The formation energy of single He 
atom at substitutional site in W is the lowest followed by that in Ta, Os and 
Re. When ZPE corrections are considered, the formation energies of He at the 
interstitial sites and substitutional sites are somewhat increased. However, 
formation energies are about 0.1 eV increased when ZPE corrections are 


included, which do not change the site-preference of interstitial He in all the 


four metals. In the following calculations, the lowest-energy configurations 
of interstitial He in the metals are taken as references. 

Table 1. The solution energy (in eV) with and without ZPE correction of 
He at various interstitial and substitutional sites in bec W and Ta, and hep 
Re and Os via first-principles calculations. The interstitial sites are plotted 


in Figure 1. 


Without ZPE With ZPE 
Ta W Re Os Ta W Re Os 
T 3.33(3.16") 6.18(6.15°) 5.55 7.19 3.41 6.26(6.14°) 5.67 7.34 
O  3.64(3.42°) 6.39(6.39") 5.29 6.64 3.73 6.43(6.37ł) 5.34 6.71 
BO — -— 4.95 6.33 -— -— 5.07 6.44 
C - -— 5.73 7.11 -— — 5.81 7.18 
Sub 1.76(1.75°) 1.62(1.38°) 2.10 1.84 1.82 1.69 2.18 1.93 


* Reference [30] 
? Reference [31] 


3.2. Diffusion of interstitial He atom in Ta, Re, and Os 

As shown in the context, the occupancy preference of an interstitial He 
atom has been determined and shown in Table 1 for Ta, W, Re, Os. The 
diffusion behavior of He from one ground-state to the next in the systems 
is determined using the nudged elastic band (NEB) method as done in our 
previous work [32, 33]. Thus, two migration paths of He atom in Ta are 
chosen, i.e., from one most-stable T site to its nearest T site, and from one 
T site to the second nearest T site passing through one O site. The barriers 
of interstitial He in perfect Ta are shown in Figure 2. We find that the most 
favorable diffusion path of one interstitial He in Ta is from T site to its closest 
T site (from T; to T2), and the barrier is as low as 0.096 eV which is in good 
agreement with the value by Yin et al [34]. The energy-lowest barrier in Ta 
is slightly larger than that in bec W and almost the same as that in bec Nb 
[35]. The barrier is 0.33 eV for He to diffuse from T site to its second nearest 
T site (from T, to T3) passing through O site. In contrast, the migration 


barrier of an interstitial He from T site to its nearest T site in perfect W is 
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even as small as 0.06 eV which agrees well with the values by Becquart [36] 
but slightly smaller than the values by González [37]. Therefore, interstitial 
He atom is relatively free to move in W than that in Ta. 


0.4 
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Figure 2: (Colour online) (a) The diffusion path of one He atom in hcp Re and Os. (b) 
The diffusion barriers of He from Tı to Tz and T3 in Ta, and barriers of He from T site 
to BO’ passing through BO, O, and C sites in Re and Os. 


In hcp Re and Os, the diffusion path and barrier are very different from 
that in bec W and Ta. As shown in Figure 2, it is difficult for one inter- 
stitial He atom to migrate from one BO site directly to another BO site. 
The diffusion paths of one interstitial He atom to move along the paths 
T—BO—>0—>C—>BO’ and T—BO—>O—BỌO’ are considered in both Re and 
Os. The diffusion barrier of one He atom along T—+BO is close to zero, and 
T site is the transition state. In Re, it has to overcome 0.31 eV for one He 
atom to move from BO site to the transition state O site, and 0.57 eV for 
He atom to diffuse from the transition state O to another transition state C 
site and then directly to BO’ site. In contrast, it will only overcome 0.31 eV 


for one interstitial He atom to move from BO to BO’ site passing through 
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transition state O site. Thus, interstitial He atom has to overcome 0.31 eV 
to move from one ground-state to another in Re. The chosen diffusion paths 
of one interstitial He atom in Os are the same as that in Re. The diffusion 
barrier of He from T to BO site is also close to zero. The diffusion barrier 
of one He atom to move along the path BO-O is 0.24 eV in Os, and it has 
to overcome 0.70 eV for the path of O+C-+BO’. However, it only needs to 
overcome 0.24 eV for He atom to diffuse from one BO to BO’ passing through 


transition state O site. 


3.8. Clustering of interstitial He in Ta, Re, and Os 


foot hd 
oo- b 


| 
Figure 3: (Colour online) The close-packed helium monolayer structure of He, cluster 


He 


(n=2-6) between (110) planes in Ta. The grey and green balls are Ta and He atoms, 
respectively. 


The segregation of He atoms is associated with cluster and bubble for- 
mation in metals, which directly influences the mechanical properties and 
performance of the metals in fusion reactor. In our previous work, it has 
been found that He atoms are energetically favorable to accumulate in a 
close-packed arrangement at interstitial sites between (110) planes forming 
helium monolayer structure in W [38]. The accumulation of multiple He 
atoms at interstitial sites in Ta, Re and Os forming He,, cluster is studied 
in the same way. Considering the symmetry of the system, all possible in- 
terstitial sites closest and next closest to any He atom in the ground-state 
He,,_; cluster are considered and the formation energies of the configurations 
are compared to find the most stable state of He„. As many as nine possible 
sites close to the tetrahedral interstitial He are calculated to find the ground- 
state configuration of Hes. The most stable He, cluster is found in the same 
way once the formation energy of the configuration is the lowest. The stable 


configuration of He, cluster in W was exhibited in our previous work [38]. 


Like W, metal Ta is also a bec structure. Single He atom prefers to occupy 
T site in Ta, and the second He atom prefers another T site, forming He-He 
pair along (100) direction as the first structure in Figure 3, which is about 
0.16 eV lower than the energy of the next stable configuration. It is found in 
recent work that He-He pair is slightly displaced in the (100) direction in bec 
Nb [35]. The third He atom will result in the displacement of the original two 
He atoms from T sites, forming a triangle between (110) planes. With the 
continuing segregation of He atoms at the interstitial sites, the close-packed 
interstitial He, cluster is formed between the (110) planes. The structures of 
He,, cluster (n=2-6) between (110) planes in Ta are shown in Figure 3. In bec 
Ta, the distance between (110) planes ay is the largest interplanar spacing, 
where do is the lattice constant of Ta. In metals Re and Os, BO is the most 
stable interstitial site for He. In hcp Re, the first He atom occupies BO site, 
and the second He atom takes another BO site along (0001) direction. The 
third He atom occupies the midway of the original two He atoms. When the 
fourth He atom diffuses to the cluster, the cluster becomes more compact. 
When the fifth He atom segregates to the He cluster, the cluster becomes 
a hexahedron, and the sixth He atom stays outside of the hexahedron. As 
for the situation of interstitial He cluster in Os, the first He atom also takes 
BO site, while the second He atom occupies another BO site in the same 
octahedron. The following He atoms do not take the regular high-symmetry 
interstitial sites but segregate together and form compact He,, cluster. 

The incremental binding energies of interstitial He, in W, Ta, Re and 
Os are calculated and plotted in Figure 4. The positive value of the binding 
energy suggests that He atoms can accumulate in the clusters one by one. The 
general tendency of the binding energy increases with the He atom number in 
all the four systems. The influence of ZPE corrections on the energies is also 
considered. It is found that the influence of ZPE corrections is very small to 
the binding energy and thus is negligible. The binding energy of He, cluster 
in W is the largest followed by that in Os, Re and Ta. The continuous increase 
of the total binding energies of He, cluster with n in metals indicates that 
trapping one more He atom to the He,, is possible. Specifically, the binding 
energies for one He atom attracted to Hes are as high as 2.35 eV and 0.85 
eV in W and Ta, respectively, where the results in W agree well with the 
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Figure 4: (Colour online) The incremental binding energies of He, clusters in Ta, W, Re 
and Os as a function of the number of He atoms in the clusters. 


value obtained by Becquart et al [39]. The binding energies of one He atom 
to Hes are 1.46 and 1.63 eV in Re and Os, respectively. So the tendency for 


He atoms to segregate in all the four metals still exists. 


3.4. Clustering of He in vacancies in Ta, Re, and Os 


In a fusion environment, the implantation of He ions and high-energy 
neutrons knocks on the lattice atoms and produces a large amount of va- 
cancy and vacancy clusters in materials. The produced vacancy and vacancy 
clusters inevitably interact with the ubiquitous He atoms, resulting in the 
formation of He clusters and bubbles, and ultimate hardening and embrit- 
tlement of the material. It is difficult for experiments to find the lattice 
location of He in W despite careful studies by Picraux et al, however, it is 
qualitatively found that multiple He atoms have the tendency to be trapped 
at defect centers [40]. In the presence of vacancy, single He atom is much 
more preferable to take vacancy center (substitutional site) rather than inter- 
stitial site, and thus vacancy can act as strong traps for He atoms diffusing 
through the bulk and form VacHe,, cluster. In this part, the energetics and 
configurations regarding the segregation of He atoms inside vacancy clusters 
are considered in W, Ta, Re and Os. The ground-state configurations of 
VacHe,, clusters are determined as that done for He». When two He atoms 


are trapped inside one Ta vacancy, the two atoms form a split dumbbell along 
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(110) direction. In contrast, two He atoms in a W vacancy reside along (111) 
direction. When the number of He atom is added up to 3 in a Ta vacancy, 
a equilateral triangle is then formed similar to the situation in a W vacancy. 
Four He atoms form a tetrahedron in monovacancy in Ta and W. When He 
atom number is further increased in Ta vacancy, He cluster becomes more 
and more compact but irregular. As for two He atoms accommodated in 
monovacancy in Re and Os, the situations are similar to that in Ta and W 
although the orientations of the two-He-atom dumbbells are very different. 
The clusters containing three and four He atoms form triangle and tetrahe- 
dron in Re and Os, respectively. The He, cluster is more compact in Os than 
that in Re. 


b m- Tae W + Re -e- Os O NAGS d g (b) 
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Figure 5: (Colour online) (a) The incremental binding energies of VacHe,, clusters as a 
function of the number of He atoms in Ta, W, Re and Os. (b) The incremental binding 
energies of Vac2He,, clusters in Ta, W, Re and Os. 


The incremental binding energies are calculated regarding the ground- 
state configurations of Vac,,He, and exhibited in Figure 5. The binding 
energetics of VacHe, are shown in Figure 5 (a). The binding energy of one 
He atom with monovacancy in Ta is 1.67 eV, and the value becomes 1.70 
eV once ZPE correction is considered. The values are slightly smaller than 
the value obtained by Yin et al [34], which may be due to the different 
convergence criteria concerning the forces on each atom. In general, the ZPE 
correction is very small and negligible for the binding energy of VacHe, in Ta. 


The similar situations are found in W, Re and Os, and ZPE corrections are 
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negligible. The binding energies of one He to vacancy and VacHe are 4.55 and 
3.21 which agree well with the results by Becquart [41] and Nguyen-Manh 
[42]. Basically, the binding energy of VacHe,, decreases with He number 
although there are some fluctuations. The calculated binding energies of 
VacHe,, in W agree well with the values by González et al [41, 43]. This 
indicates that the binding ability of He atom to VacHe,, is reduced gradually 
with n. The underlying physical reasons controlling the binding energetics 
of He with vacancy will be analyzed in the following part. In the same way, 
the binding energies of He atoms with two-vacancy cluster (Vac2He,,) in Ta, 
W, Re and Os are calculated and presented in Figure 5 (b). The similar 
decrease tendency is found for the binding energy of VacjHe,,. According to 
the present calculations, the binding energy of Vac,,,He, is much larger in W 
than that in Ta and Re especially in Ta. The binding strength of VacmHen 
in W is close to that in Os. 


3.5. Mutation of He cluster in Ta and W 


In order to investigate the growth and evolution of a He, cluster, it is 
necessary to study the possibility for the emission of lattice atoms around 
the interstitial He, in Ta, W, Re and Os systems. The present calcula- 
tions suggest that the self-interstitial atom in (111) dumbbell is the most 
stable configuration in both W and Ta. The ground-state configurations 
of SIA in Ta and W agree well with the results obtained by Nguyen-Manh 
et al [44]. Therefore, it is assumed that lattice atoms are emitted from 
the He, cluster and the lattice atoms move away by (111) dumbbell con- 
figuration in Ta and W, respectively. The He, cluster can be dissolved by 
the reactions of He,—VacHe,+SIA or He,—He,_;+He, where SIA suggests 
self-interstitial atom. The reaction of He,—->VacHe,,+SIA can be defined by 
Eq = E(VacHe,,) + Esgra— E(He,) — E, and the reaction of He,—He,_1+He 
can be defined as Eq. 2. The dissolution energy can defined by E)=Eqt+Em, 
where E,, is the diffusion barrier of interstitial He atom or SIA in perfect 
system [22]. The defined dissolution energy suggests the energy needed for a 
SIA or He atom to emit from their initial site and diffuse away. If the disso- 
lution energy of a SIA from the He,, is smaller than that of a He atom from 


the cluster, it will be possible for the mutation of He, to VacHe,. Otherwise, 
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it can only emit He atom from He,. According to our results in the context, 
the diffusion barriers E,,, for an interstitial He atom from one T site directly 
to its nearest T site are 0.06 eV and 0.096 eV in W and Ta systems, respec- 
tively. The diffusion barriers of a SIA translated from one (111) dumbbell to 
the nearest (111) are 0.003 eV and 0.027 eV in Ta and W, respectively. 
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Figure 6: (Colour online) The dissolution energies of SIA and He atom from He,, cluster 
as a function of He atom number in Ta and W. 


In bec Ta and W, the dissolution energies of SIA and He atom from Hen 
are calculated and plotted in Figure 6. From the Figure, the dissolution 
energies of SIA from He, cluster in Ta and W are initially larger than that 
of He atom from He, clusters. The dissolution energies of a SIA from He, 
decrease when n is increased in both W and Ta. In contrast, the dissolution 
energies of a He atom from He, cluster increase with n in both systems. 
When the number of He atom n exceeds 6, the dissolution energies of a He 
atom from He, cluster become larger than that of a SIA from He, in both 
Ta and W. The results indicate that a SIA can be emitted from He,,, forming 
VacHe,, cluster in both systems. Similarly, in both Re and Os, the dissolution 
energies of a SIA from He, cluster decrease with n, while the dissolution 
energies of a He atom from He, increase with n. From the calculations, it 
has been demonstrated that the emission of a SIA from He,, to form VacHe,, 
and Vac,,He, is possible if n is large enough. By the continual emission 
of SIA and trapping He atoms, He, cluster may grow larger and form He 


bubbles. However, there is no systematical difference between Ta and W for 
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the grow ability of Hen. 


4. Discussion 


4.1. The physical reasons controlling He clustering in Ta, W, Re, and Os 


The underlying mechanisms controlling the clustering of He atoms in Ta, 
W, Re and Os are analyzed in this section. As is known, He is a kind of 
closed-shell atom, and any hybridization of He with Ta, W, Re and Os is 
energetically unfavorable. Therefore, He atom prefers to accumulate at free 
space to avoid the hybridization with bulk atoms in the systems. He atom at 
interstitial site often results in the displacement of surrounding lattice atoms 
away from the He atom. For example, He atom located at T site in Ta pushes 
the closest lattice atom away by 0.19 A. He atoms at BO sites in Re and Os 
lead to the displacements of the nearest lattice-point atoms away by 0.24 and 
0.21 A, respectively. In this way, He can stay away from the lattice atom to 
avoid hybridization. This may be the reason why the binding energy of He,, 
increase with He atom number. In contrast, vacancy and vacancy cluster can 
provide more free space for He atoms to avoid chemical hybridization with 
metal atom. As the increase of He atom in vacancy cluster, the free space 
in the vacancy cluster is gradually reduced which result in decrease of the 
binding energy of He with Vac,,Hen. 

The volume change of the system may also associate with the incremen- 
tal binding strength of He clusters. The volume change of the system con- 
taining a He, is defined as AV=V(He,,)-V(Perf), and the volume change 
of the system containing a VacHe,, is approximated by AV=V(VacHe,,)- 
(V(Perf) /128)*127. The relationships regarding the binding energies of He, 
and VacHe,, clusters with the volume change of the systems containing the 
clusters are calculated and presented in Figure 7. The AV of the system 
containing He, cluster increases with n. From Figure 7 (a), the incremental 
binding energies of He, increase with AV of the system containing the cluster 
in all the four metals. The situation of VacHe,, is different from that of Hen. 
He atoms are energetically favorable to cluster close to vacancy center. If the 
space of vacancy is sufficient large, and He atoms in vacancy may have less 


influence on surrounding lattice atoms. The AV of the systems containing 
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Figure 7: (Colour online) The relationship between the volume change and the incremental 
binding energy of Hen (a) and VacHe,, cluster (b). 


VacHe,, clusters also increases with n, while the binding energies decrease 
with AV. Although the accumulation of He atoms inside vacancy can result 
in volume increase, the increased space of single vacancy still can not meet 
the space-requirement of He atom. 

The accumulation of He atoms at interstitial sites forming He, induces 
the increase of free volume in the systems, which also is also associated 
with lattice distortion around He,, cluster. The lattice distortion induces 
local stress, which reduces the vacancy formation energy around He,. The 
decrease of vacancy formation energy suggests the stability of lattice-point 
atoms are substantially reduced. From the definitions of vacancy formation 
and dissolution energies close to He, cluster, the difference is only a constant. 


This is why both vacancy formation and dissolution energies are reduced with 
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He atom number n. The reduction of the stability of lattice-point atom may 
be the underlying reasons why the emission of a SIA from He,, clusters is 


possible. 


4.2. Comparison regarding the nucleation and evolution of He clusters in Ta, 
W, Re, and Os 

In a fusion reactor, neutron irradiation results in transmutation of W to 
Ta, Re and Os isotopes [2], and the transmutation elements segregate in W 
forming precipitates containing high concentration Ta, Re and Os [12, 17- 
19]. As the continual segregation of transmutation elements, W alloy may 
transform to bcc+o, g, x and eventually bec Ta, and hcp Re and Os systems. 
The precipitates naturally interact with He and influence He retention and 
He bubble formation in W. The interaction of He plasma with W often re- 
sults in bubble and blister formation [45-47], which degrades the mechanical 
properties of material and threatens the stability and safety of reactor. It 
has been widely observed that the presence of transmutation elements Ta, 
Re and even their precipitates can reduce D retention and suppress bubble 
formation [8-11]. However, the influence of transmutation elements Ta, Re 
and Os on He retention and bubble formation involving the behaviors of He 


in pure Ta, Re and Os systems is mysterious. 
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Figure 8: (Colour online) The stability of isolated He, clusters (n=2 to 6) in Ta and W 
calculated using objective kinetic Monte Carlo method in the same conditions. 


According to the present calculations, He atoms are energetically favor- 


able to segregate in Ta, Re and Os similar to the situation in W. The differ- 
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ence is that He can continually accumulate at interstitial sites in W and Ta 
forming closed-packed monolayer between (110) planes, whereas He atoms 
form disordered structures in Re and Os. The binding strength of He, clus- 
ter in W is the strongest followed by that in Os, Re and Ta. Based on 
the energetics and kinetics properties of He in W and Ta, the desorption of 
isolated interstitial He clusters in the systems is calculated using objective 
kinetic Monte Carlo method and presented in Figure 8. The heating rates 
is 1 K/s. It is found that the desorption peaks of interstitial He, clusters 
mainly reside between about 140 and 340 K in Ta. In contrast, the desorp- 
tion peaks of He, clusters are shifted to the range of 410 to 850 K in W. So 
He,, clusters are much more stable in W than in Ta. In a word, the larger He 
clusters is more stable due to larger binding energy. The presence of stable 
Hen clusters causes the emission of self-interstitial atom from the clusters, 
that is, trap-mutation process occurs. The trap-mutation process goes on 
and may result in bubble formation. The migration of He atom in W and Ta 
is easy due to the small diffusion barriers. In contrast, the diffusion of single 
He atom in Re and Os is relatively difficult because the diffusion barriers for 
the favorite path in Re and Os are larger than that in Ta and W. Generally, 
all these conditions may contribute to the suppression of Ta, Re and Os on 
He bubble nucleation and growth in W. 


5. Conclusions 


Systematic first-principles calculations are performed to study the occupy, 
diffusion, and segregation of He atom, as well as the underlying evolution 
mechanism regarding He clusters in bcc Ta and W, and hcp Re and Os. 
It is found that He atom prefers tetrahedral (T) interstitial site in Ta, and 
basal octahedral (BO) interstitial site in Re and Os. The diffusion study of 
interstitial He using NEB method shows that the favorite diffusion paths of 
one He atom are from one tetrahedral (T) site to the closest T in Ta. The 
barrier for the favorite path is 0.096 eV in Ta. He atom is more preferable 
to move from one basal octahedral (BO) site to the next passing through 
octahedral site along (0001) direction in both Re and Os, and the barriers 
for the diffusion in Re and Os are 0.31 and 0.24 eV, respectively. He atoms are 


energetically favorable to segregate at interstitial sites forming He,, cluster in 
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all Ta, W, Re and Os, and He atoms are energetically favorable to segregate 
at interstitial sites forming close-packed monolayer structure between (110) 
planes in bcc Ta. The presence of vacancy promote the segregation of He 
atoms in Ta, W, Re and Os. He atoms aggregate in vacancy cluster forming 
Vac He, in the systems. The binding strength of He, and VacmHen clusters 
in W is relatively large than that in Ta, Re and Os. The presence of solutes 
Ta, Re and Os, especially Ta, may contribute to suppress the He bubble 


nucleation and growth in W. 
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